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Введение. Рассматривается метод решения задачи о действии вертикальной силы, движущейся по 
произвольному закону, на бесконечную пластину. Данный метод и полученные результаты могут 
использоваться при изучении действия подвижной нагрузки на различные конструкции. 

Материалы и методы. Разработан оригинальный метод решения задач о действии на бесконечную пластину, 
лежащую на упругом основании, вертикальной силы, движущейся произвольно по незамкнутой кривой 
произвольной формы. Для этого используется фундаментальное решение дифференциального уравнения 
динамики пластины, лежащей на упругом основании. Считается, что движение силы начинается в достаточно 
удаленный момент времени. Поэтому начальные условия в такой постановке задачи отсутствуют. При 
определении фундаментального решения выполняется преобразование Фурье по времени. При обращении 
преобразования Фурье изображение раскладывается по параметру преобразования в ряд по полиномам Эрмита. 
Результаты исследования. Представлено решение задачи о бесконечной пластине, лежащей на упругом 
основании, по которой с переменной скоростью движется сосредоточенная сила. В качестве траектории 
рассматривалась гладкая незамкнутая кривая, состоящая из прямых и дуг окружностей. Изучается поведение 
компонент вектора перемещений и тензора напряжений в месте расположения движущейся силы, а также 
процесс распространения энергии волн, для чего рассматривается изменение вектора плотности потока энергии 
Умова-Пойнтинга. Исследуется влияние на перемещения, напряжения и распространение упругих волн 
скорости и ускорения перемещения силы. Изучается влияние формы траектории движения силы на 
напряженно-деформированное состояние пластины и на характер распространения упругих волн. Результаты 
свидетельствуют о том, что метод достаточно устойчив в широких пределах изменения скорости движения 
силы. 

Обсуждение и заключения. Проведенные расчеты показали, что наиболее существенным фактором, влияющим 
на напряженно-деформированные состояния пластины и на распространение энергии упругих волн вблизи 
сосредоточенной силы, является скорость ее движения. Данные результаты будут полезны при изучении 


динамических процессов, порождаемых подвижной нагрузкой. 
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Введение. Закономерности динамических процессов в твердых средах, вызванных действием 
подвижной нагрузки, представляют значительный интерес, а решения подобных задач находят многочисленные 
приложения и предполагают использование разнообразных методов. В ряде работ для исключения времени из 
числа независимых переменных вводилась подвижная система координат [1-2] или рассматривалась 
квазистатическая постановка задачи [3—6]. При решении указанных задач достаточно эффективными оказались 
метод конечных элементов [7], вариационные [8—10], а также прямые методы [11-13]. В работах [14—15] 
применялся метод граничных интегральных уравнений, а в работе [16] — метод, основанный на применении 
фундаментальных решений соответствующих дифференциальных уравнений. В настоящей работе этот метод 
применяется для решения задачи о действии на бесконечную пластину, лежащую на упругом основании, 
вертикальной силы, движущейся по незамкнутой кривой произвольной формы. 

Постановка задачи. Следуя [17, 18], указанная задача сводится к решению уравнения: 


а Р 
АИ +с9 И +КИ =—, (1) 
р) 
ЕНЗ 
где О — прогиб пластины; О = ео Е — модуль Юнга; и — коэффициент Пуассона; Н — толщина 
и 
9 Н К, 
пластины; с“ =7>; р — плотность материала пластины; К = . К, — коэффициент жесткости упругого 
основания. 


Решение данного уравнения соответствует потоку энергии, направленному от источников возбуждения 
в бесконечность. Будем полагать: 


Р = 6(х - хо(Е))6 (у — Уо (0). 


Данная сила перемещается по незамкнутой траектории \, начало и конец которой уходят в бесконечность. 
= х(:) 


У=у(!) 


Параметрическое задание траектории имеет вид: | ‚ где { — время. Считается, что сила начинает 
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движение в начале траектории, расположенном на достаточном удалении от того места, в котором исследуется 
ее воздействие на пластину в момент времени = —ю. Поэтому начальные условия в такой постановке 
отсутствуют. 


Материалы и методы. Рассмотрим фундаментальное решение уравнения (1), которое можно получить 
из уравнения: 


ДРИЙ + с-20РИ/ + КИ! = 5(х — хо) (у — Уь (Е — т. (2) 
Известно, что решение уравнения (1) можно представить в виде: 
бк, в = [ДИ хо, У У — ОР(х у, ахоауь ах. 
В нашем случае с учетом конкретного вида движущейся силы имеем: 
ИЖ У,В) = Г, И (, хо (®), У, Уо (Е — тат. 
Применяя преобразование Фурье по времени к уравнению (2), получим дифференциальное уравнение: 
ДРИ, — РС, + ИБ = 56(х — хо) 6 — узде". (3) 
Используя принцип предельного поглощения и преобразование Фурье по переменным х и у, и при 


2 
[99 
условии К > 2 Можно получить решение уравнения (3): 


2 г. 
Ио (хо, У) = тео [Ко(ваВ) — Ко(а» В), 


где А = (хх у +(5-у, у" = || _ м - аа = Хе" /4; а» = Хей /4; К, (2) — функция Макдональда. 


02 
При условии К < > решение уравнения (3) выглядит следующим образом: 


"НО (ХВ) — К В)|, 


Ио (х, хо, У, Ус, 6 65? /с^) = ы 


2 
где х = К —К; НГ” (2) — функция Ганкеля. 


Для обращения преобразования Фурье решение И (х, хо, у, Уд, о? /с?) раскладывается по переменной 


4 т 


2 
(0) в _@ (69) 
— в ряд по системе ортогональных функций | Ис? Н К (©) где Н,(2) — полиномы Эрмита. 
С 


2 
(69) в (0 
Учитывая, что функция И (х Хо, У, У») является четнои по р ‚ В разложении будут присутствовать 


только четные слагаемые. Тогда: 
2 
И 2 с?) = у сан (° 
о(х, Хо, У, Уо, 9 /С^) = Ук-ом”2к (Х, Хо, У, Уо)е с П?2к р ‚› где 


1 со ый. 
мак (Х, Хо, У, Уо) = Брак 5 Мо, Хо, У, У 2”) е "Нок (рае. 


Учитывая соотношение: 


со 2 | С2Е2 
| г" Ис? Ноу (©) е ‘и = 2с | (-1)Ке- "2 Нк (се), 


получим: 


22 _ С2Е2 
Ио (х, Хо, У, Ус, Е) - 2 Уи быжыь уже р /2Нок(сь). 


В этом случае решение исходного дифференциального уравнения будет иметь вид: 


п со со К _с*@-9"/ 
Их у, в) = 2с |^2 [® У (Тик (хо (©, у, 5 6е 2Нок(с(Е — ©) а. 
Сделав замену переменной интегрирования, получим: 


В $2 $2 :)242(- 1 
Обл) => | | ох кА ат СТАТЫ 


Такой вид решения позволяет для вычисления интеграла применить квадратурную формулу Гаусса-Эрмита. 

Чтобы улучшить сходимость ряда использовался метод Куммера. Следуя этому методу, необходимо 
подобрать ряд, сумма которого известна, а разность исходного ряда и ряда подобранного должна представлять 
быстро сходящийся ряд. В качестве такого ряда можно взять: 


ув = У [И се, ж(е - 52), у.е 52/5), 4) 
где Я — некоторая неотрицательная величина. 
Проинтегрировав по переменной т и просуммировав, получим: 


0” (х, у, Е) — псу2 И (х, Хо (Е), У, Уо (6), 4). 
Окончательно для решения уравнения (1) получим следующее выражение: 


ЭН, (82. )Н ‚(12 дат. 


2/2(-—=1)^ 
(2к)!22К 


е` "+ Нок (5/2) Нок (<У/2)а5ат, 
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И(х, ув) = Ц*(х, ув) + 


+) || [обежо@ — УИ оуруь@ — Уолт)  иобьжь(е — У /у,уз(е  5У2/5,9))х 


К=о 


2/2(—1)^ 
х ес ни (5 )ны(тУВ)азвх 
Для суммирования ряда применялся способ средних арифметических. При этом полагалось: а = 0. 

Определив прогибы пластины, можно по известным формулам вычислить остальные компоненты 
вектора перемещения и тензора напряжений в любой ее точке. Для анализа перемещения энергии упругих волн 
в пластине вычислялся вектор плотности потока энергии Умова- Пойнтинга: 

Е = —(6,й + в, -— (вхуй + 6,5). 
Результаты исследования. Проведены расчеты для случая, когда сила двигается по траектории, 


состоящей из прямых и дуг окружностей (рис. 1). При этом принимались следующие значения параметров: 
2; -4 

Н= 0,25 м; с=221 м/с; Е= 232469 Н/м” и=0,36; К= 1,864 м’. Параметры закона движения силы по 

траектории подбирались таким образом, чтобы в рассматриваемый момент времени сила всегда находилась в 


одной и той же точке траектории, отмеченной звездочкой, имея различные значения скорости у и ускорения а, а 


также при разных значениях радиуса траектории В›. Для изучения напряженно-деформированного состояния 
пластины вычислялись перемещения и напряжения вблизи точки приложения силы. 


У, м 
15 


10 





-20 -15 -10 -5 0 е 10 15 Х, м 


Рис. 1. Траектория движения сосредоточенной силы 


На рис. 2, 3 представлено изменение перемещений и напряжений при движении сосредоточенной силы 


ь 2 
по заданной траектории при у= 25 м/с, а=0 м/с”, К. =5м. Изменение указанных величин вдоль оси У 
практически не отличается от их изменения вдоль оси Х. 











У, У, О, мх10” 
Г на 2 ше т ва Ра т — ег | —] 
. 74 \М/ 
Ф-® | С 
| Е 
-0,5 Г / 4 
= \ / 
| | 
-1,0 \ \ | | 
\ Г 
я -1,5 | \ / | 
\ ] 
= 4 
[ев -2.0 г \ 1 1 
> = } \/ ‹ 
< 2,5 \ 
с 
Е -3,0 м к ЗИ», ЗИ и, ЗИ р 98 ИЕ А. 2 





— 


-8 -6 -4 -2 0 2 - 6 Х, м 


Рис. 2. Изменение перемещений: \! — вертикальное; Ц — вдоль оси Х; У — вдоль оси У 
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Зу,Зх,5ху, Н/м” 





Рис. 3. Изменение напряжений: 
\/ — вертикальное; Ч — вдоль оси Х; У — вдоль оси У 


На рис. 4 представлено перемещение энергии упругих волн вблизи сосредоточенной силы, положение 
которой на траектории обозначено красной точкой. Векторы определяют количество и направление переноса 


энергии в данной точке. 





Х, м 
Рис. 4. Вектор плотности потока энергии при у = 25 м/с, а = 0 м/с”, В =5м 
Расчеты показали, что при увеличении скорости движения силы качественного изменения 
перемещений и напряжений не происходит, а происходит лишь их количественный рост. Незначительное 
изменение качественного поведения перемещений и напряжений наблюдается только при достаточно больших 


о 
скоростях, когда выполняется условие у > с. Это следует из рис. 5, 6 (а = 0 м/с", В =5м). 
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У, мх1 07 


2 





Рис. 5. Изменение вертикальных перемещений при у = 275 м/с 


У, мх1 07 





1 0 ниваРА- 1 0 


У, м -20 
Рис. 6. Изменение вертикальных перемещений при у = 75 м/с 


На рис. 7, 8 представлено изменение максимальных вертикальных перемещений У и напряжений 5, , 5, в 
зависимости от скорости движения силы при а=0м/с’, В. =5м. Остальные компоненты перемещений и 
напряжений принимали достаточно малые значения и поэтому не представляли конструктивного интереса при 
анализе напряженно-деформированного состояния пластины. 

Вычисления, проведенные при различных значениях ускорения и радиуса В›, показали, что эти факторы 
слабо влияют на напряженно-деформированное состояние пластины. Также слабо зависит от указанных факторов 
качественная картина распространения энергии волн вблизи сосредоточенной силы. 
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244 Рис. 7. Зависимость максимальных вертикальных перемещений от скорости движения сосредоточенной силы 
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Рис. 8. Зависимость максимальных напряжений от скорости движения сосредоточенной силы 


Обсуждение и заключения. Наиболее существенное влияние на напряженно-деформированное 
состояние пластины и распространение энергии упругих волн вблизи сосредоточенной силы оказывает 
скорость ее движения. Радиус кривизны траектории и ускорение движения силы влияют несущественно. 

Результаты расчета свидетельствуют о том, что метод решения задач о действии подвижной нагрузки 
достаточно устойчив в широких пределах изменения скорости ее движения. Метод экономичен и прост, т. к. в 
нем используются уже известные фундаментальные решения. 


Библиографический список 

1. Опаг баба. Кезропзе оРа ЗО @азис Ба -зрасе © а 41516 \{е тоушо 1оа4 / Опиг За, Ваг1$ Егфа$, 
Ма! Езе // Насейере Логпа| ог Мафетанс$ ап 5‘ай$с$. — 2017. — Уч. 46 (5). — Р. 817-828. 
В рз://Ао1.0го/10.15672%20/Н]М$.2017.434 

2. Динамика слоистого полупространства под действием движущейся и осциллирующей нагрузки / 
В. В. Калинчук, Т. И. Белянкова, Г. Шмид, А. Тосецки // Вестник Южного научного центра РАН. — 2005. — 
Т. 1, №1. — С. 3—1. 

3. Уцуап СВеп. Оупапис гезропзе оЁ ап еазяс рае оп а сгоз$$-ап1зогор1с рогоеЙазис Ва Мр1апе №0 а 10а 
тоуе оп 15 за асе / Упуап СВеп, МПа О. ВезКои, Лап® Олап // ЗоЙ Оупапис$ апа Еапачаке Епошеети®. — 
2018. — Уа1. 107. — Р. 292-302. В@рз://401.0г®/10.1016/.зоПауп.2018.01.038 

4. Карапоу, УТ. ТБе еее \ауе оп ап @азисаПу зирроцед КисЬвоЁ ре / лаз Каршпоу, РапПа 
А. Рика7тсКоу, ОтаВат А. Козегзоп // ТВе Лоигпа! оЁ фе Асочзйса| Зослеу ог Атепса. — 2014. — Уо]1. 136 (4). 
— Р. 1487—1490. В@рз://4о1.оге/10.1121/1.4894795 

5. Егорычев, О. О. Воздействие подвижной нагрузки на многослойную вязкоупругую пластину, 
лежащую на вязкоупругом основании / О. 0. Егорычев // Вестник Московского государственного 
строительного университета. — 2007. — Вып. 1. — С. 39—42. 

6. Динамическое поведение безграничной упругой пластинки при воздействии подвижной (бегущей) 
нагрузки / М. Ж. Досжанов, Е. Н. Искак, Б. Ж. Сактаганов [и др.]| // Путь науки. — 2016. — Т. 1, № 11 (33). — 
С. 26—28. 

7. Шишмарев, К. А. Постановка задачи о вязкоупругих колебаниях ледовой пластины в канале в 
результате движения нагрузки / К. А. Шишмарев // Известия Алтайского государственного университета. — 
2015. — № 1/2 (85). — С. 189-194. Вирз://Ао1.оге/10.14258/17уази(2015) 1.235 

8. Пушемлтст, В. Угайоп$ оР а Мш@аПа р] ме заб]е‹е ю а раш оЁ шегйа| 10а4$ шоуше ш оррозие 
фтесяоп$ / Ваг4опте] Рутеум1с7т, ошииК Разаг$К1, Схеза\ Г. Ва]ег // оигпа[ о? Зоипа апа Утгайоп. — 2017. — 
\Уо1. 386. — Р. 265—282. В@рз://401.0ге/10.1016/.]зу.2016.09.027 

9. Езеп, Г. А пем Нице е@етепе Гог гапзуегзе уганоп оЁ тесбапошаг @ш рез ипдег а тоуше таз$ / 
15та! Езеп // Ешце Еетепз ш Апаууз$ апа Оезюп. — 2013. — Уча. 66. — Р. 26-35. 
6 рз://401.ого/10.1016/.йпе1.2012.11.005 

10. оп Ошевиа. У1гаНоп апа[у$1$ оР иснопаПу отаде4 р]=ме ут а тоуше та$$ / Зопе Ошеваа, Ш 
ЛаВао, ли 7Ваполап® // АррПеа Мафетайса! Мо4деШши». — 2017. — Уч. 46. — Р. 141-160. 


6 р5://401.0г2/10.1016/.арта.2017.01.073 


Механика 


245 


БИр://уезК-допза.га 


246 


А4уапсеа Епятеетиз Кезеагсй 2021. Т. 21, №3. С. 239—246. 155М 2687—1653 





11. Рагатейлс ау ог дупапис гезропзе оЁ зап4уясВ ре ипдег тоуте 10а4$ / Ошевиа Зоп2, /ваполап® 
ли, МЛабао Зш, У! \Уап // Тып-\УМаШеа Зиоасигез5. — 2018. — Уа. 123. — Р. 82-99. 
В рз://Чо1.ого/10.1016/.6%5.2017.11.012 

12. Тите-дотала $гасбига|-асоизИс апа[у$1$ оР сотрозце р|аж(ез заб ]ее © тоуше 4упаплс 1юа4$ / Уегао 
Оц, У\Уепшше ИВапе, Г/Бщке Реп, Счапе Мепо // Сотшрозие ЗЭгасагез. — 2019. — Уо1. 208. — Р. 574—584. 
6 рз://401.0г2/10.1016/.сотпрзгис{.2018.09.103 

13. Еоуоч7а%, М. А. Ап апауйса|-питегса1 зо[айоп © а5зез5 фе упатс гезропзе оЁ у15сое|азис р|а{ез о а 
тоу1е та$$ / М. А. Еоуоц7а%, Н. Е. ЕеКкапс 1, М. МойА // АррПед Мафетайса1 МодеШие. — 2018. — Уа1. 54. 
— Р. 670—696. №рз://401.0г®/10.1016/.артт.2017.07.037 

14. Галабурдин, А. В. Применение метода граничных интегральных уравнений к решению связных 
задач термоупругости с подвижной нагрузкой / А. В. Галабурдин // Известия вузов. Северо-Кавказский регион. 
Естественные науки. — 2012. — №4. — С. 29-31. 

15. Галабурдин, А. В. Применение метода граничных интегральных уравнений к решению задач о 
движущейся нагрузке / А. В. Галабурдин // Известия вузов. Северо-Кавказский регион. Естественные науки. — 
2015. — № 1. — С. 9—11. 

16. Галабурдин, А. В. Задача о бесконечной пластине, нагруженной нормальной силой, движущейся по 
сложной траектории / А. В. Галабурдин // Вестник Донского государственного технического университета. — 
2019. — Т. 19, №3. — С. 208—213. В@рз://4о1.оге/10.23947/1992-5980-2019-19-3-208-213 

17. Бабаков, И. М. Теория колебаний / И. М. Бабаков. — Москва : Наука, 1968. — 560 с. 

18. Бреббия, К. Методы граничных уравнений / К. Бреббия, Ж. Теллес, Л. Вроубел. — Москва : Мир, 
1987. — 524 с. 


Поступила в редакцию 26.07.2021 
Поступила после рецензирования 09.08.2021 
Принята к публикации 09.08.2021 


Об авторе: 

Галабурдин Александр Васильевич, доцент кафедры «Математика и информатика», ФГБОУ ВО 
«Донской государственный технический университет», (344003, РФ, г. Ростов-на-Дону, пл. Гагарина, 1), 
кандидат физико-математических наук, доцент, ОКСО: Бр://огс14.0г2/0000-0003-0411-6724, 


Саабигаш(@лта!.га 


Автор прочитал и одобрил окончательный вариант рукописи. 


